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ABSTRACT 

We estimate the dust temperatures of the clumps in the p Oph main cloud taking 
into account the 3D geometry of the region, and external heating from the interstellar 
radiation field and from HD147879, a nearby luminous B2V star, which is believed to 
dominate the radiation field in the region. We find that the regions where prestellar 
cores are observed (i.e. at optical visual extinctions > 7 mag) are colder than ~ 
10—11 K. These dust temperatures are smaller than those which previous studies 
of the same region have assumed. We use the new dust temperatures to estimate the 
masses of the prestellar cores in the p Oph main cloud from mm observations, and 
we find core masses that are larger than previous estimates by a factor of ^ 2 — 3. 
This affects the core mass function (CMF) of the region; we find that the mass at 
which the core mass spectrum steepens from a slope a ~ 1.5 to a slope a ~ 2.5 has 
moved from ^ 0.5 Mq to 1 Mq. In contrast with the CMF in other star forming 
regions (e.g. Orion), there is no indication for a turnover down to the completeness 
hmit (- 0.2 Mq), but the CMF may flatten at around 0.4 Mq. 

We generalize our results to the prestellar cores in Taurus and in Orion. In Taurus 
the ambient radiation field heating the prestellar cores is believed to be weaker than 
than that in p Oph. Hence, the dust temperatures of the cores in Taurus are expected 
to be below ~ 10—11 K. In Orion the radiation field is believed to be 10"^ times stronger 
than the standard interstellar radiation field. Based on this assumption we estimate 
that the dust temperatures of the prestellar cores in Orion are around ^ 20 — 30 K. 

Key words: Stars: formation - ISM: clouds-structure-dust - Methods: numerical - 
Radiative transfer 



1 INTRODUCTION 

Prestellar cores are condensations in molecular clouds that 
are on the verge of collapse or already collapsing (e.g. Myers 
& Benson 1983; Ward- Thompson et al. 1994, 1999; Ward- 
Thompson, Andre & Kirk 2002). They represent the first 
phase in an evolutionary model of star formation that is 
based on observations of different types of young objects: 
starless core/prestellar core — » Class Class I Class II 
Class III (Lada 1987; Andre et al. 1993; Andre et al. 2000; 
Di Francesco et al. 2007; Ward- Thompson et al. 2007). The 
study of starless/prestellar cores is important in constrain- 
ing the initial conditions for star formation. Prestellar cores 
have been observed both in groups (as in e.g Ophiuchus, 
Taurus, Perseus) and in relative isolation (e.g. B68; Alves et 
al. 2001). They are very cold; their temperatures are below 
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20 K and most probably around 10 K (Evans et al. 2001; 
Stamatellos & Whitworth 2003a,b; Stamatellos et al. 2004). 
They are also very dense; their typical central densities are 
^ lO^cm"^. They are observed either at NIR wavelengths 
where they are seen in absorption against the luminous back- 
ground (e.g. Bacmann et al. 2000), or at FIR and submm 
wavelengths, where they appear in emission (e.g. Motte et 
al. 1998; Ward-Thompson et al. 2002; Kirk et al. 2005). The 
peak of their emission is around 150-200 |i.m, consistent with 
the fact that they are very cold. 

Submm and mm observations are often used to deter- 
mine the masses of cores (e.g. Motte et al. 1998; Nutter et 
al. 2007). At these wavelengths the core is optically thin to 
the radiation it emits, hence the observed flux from the core 
is F\ = Tx B\{TduBt)- Therefore the column density N{H2) 
along the line of sight is 



prUH An KA -BA(rdust) 
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where Afi is the solid angle of the telescope beam for a re- 
solved source, or the solid angle of the source if unresolved, 
N(H2) is the hydrogen column density, k\ is the dust opac- 
ity per unit mass, and Tdust the temperature of the dust 
(Hildebrand 1983). The above relation assumes that the 
dust is isothermal along the line of sight within the core. 
The mass of the core is then determined from the column 
density, Mcoro ~ J N{H2)dS , where the integral is over the 
projected area of the core. The main uncertainties in de- 
termining core masses using the above method come from 
our limited knowledge of the properties of the dust in cores 
(kx) and the dust temperature (Tdust). Eq. [T] is very sen- 
sitive to the temperature, since at these low temperatures 
the Planck function is non-linear and the Rayleigh- Jeans 
approximation is not valid. Hence, even underestimating or 
overestimating the core temperature by a few degrees may 
lead one to overestimate or underestimate the Planck func- 
tion (and consequently the core mass) by a factor of 2 to 3 
(Stamatellos & Whitworth 2005b). 

It is important then to investigate whether the core tem- 
peratures estimated by previous authors using the core SED 
(defined by observations at only a few wavelengths) are ac- 
tually representative of the dust temperatures in cores. The 
aim of this paper is to use detailed 3D radiative transfer 
modelling, taking into account the core environment (e.g. 
nearby luminous sources, ambient cloud), to estimate the 
temperatures in prestellar cores. We will focus our study on 
the cores in the p Oph main cloud, but we will attempt to 
generalise our results to other star forming regions. 



2 THE p OPH MAIN CLOUD 

p Ophiuchi is a star forming region where many prestellar 
cores have been observed (see Motte et al. 1998; Nutter et 
al. 2006) along with more evolved protostars (Class 0, I, II 
objects; e.g. Motte et al. 1998; Bontemps et al. 2001; Wilking 
et al. 1989; Andre & Montmerle 1994). It is one of the closest 
star-forming regions, being at a distance from 140 to 160 pc 
(Bontemps et al. 2001; Motte et al. 1998). 

We shall confine our study to the p Oph main cloud 
(i.e. L1688), where six major clump^ have been identified 
(Oph-A, Oph-B, Oph-C ,Oph-D, Oph-E and Oph-F). The 
prototypical Class object VLA1623 (Andre et al. 1993) 
is located in the Oph-A clump. Each of these clumps con- 
tains a few tens of solar masses and has an extent of ~ 0.3 
pc. These clumps show substructure; a large number (jj 60) 
prestellar cores have been identified in them. The core mass 
function in the p Oph main cloud is similar to the stellar 
IMF, which suggests that the masses of stars are determined 
by fragmentation at a very early stage (Motte et al. 1998). 

The external radiation field incident on these clumps is 
believed to be dominated by HD147889, a nearby B2V star, 
which is located at a distance 0.5 — 1 pc away from from the 
far side of the Oph-A clump (Liseau et al. 1999). This is a 
luminous star of 5500 Lq (Wilking et al. 1989). Assuming a 
distance of 1 pc and a clump radius of 0.2 pc, if there is no 

^ To avoid confusion wc note that we use the term clump to refer 
to the larger structures and the term core to refer to the conden- 
sations, i.e. the substructure, in these clumps. Some authors (e.g. 
Motte et al. 1998) use the reverse terminology. 



attenuation, 1% of the star's radiation will heat the clump, 
i.e. 55 Lq. This radiation dwarfs any radiation from other 
sources (e.g. the interstellar radiation field, nearby young 
protostars). 

Observations of the wider area of the main p Oph cloud 
(~ 4 deg^) using SCUBA on the JCMT and comparison 
with 2MASS extinction maps as part of the COMPLETE 
project (Ridge et al. 2006) reveal that there are no prestel- 
lar cores (at least down to the SCUBA detection limit) at 
visual extinctions < 7 — 10 mag (Johnstone et al. 2004). It 
is presumed that the UV radiation present in regions of low 
visual extinction makes these regions hostile for star forma- 
tion (McKee 1989). Surveys in other star forming regions 
have also revealed similar extinction thresholds that appear 
to be connected to the environment of the region (e.g. in 
Perseus the extinction threshold is lower. Ay > 5 — 7 mag; 
Kirk et al. 2006). 

Thus, it appears that most cores are quite embedded 
in their parent molecular clouds. Previous radiative transfer 
simulations of embedded cores suggest that such cores are 
very cold (Stamatellos & Whitworth 2003a, b; Stamatellos 
et al. 2004). The temperatures at their centres is around 
6 K and at their edges is from around 11 K (for a core em- 
bedded in a parent cloud of Ay = 5) to 9 K (for a core 
embedded in a parent cloud of Ay — 20). This is because 
the deeper the core is embedded the lower is the strength 
of the radiation field that heats the core. Hence, it is ap- 
propriate to ask whether previous studies of prestellar cores 
have assumed core temperatures that are consistent with the 
fact that these cores are deeply embedded in their parental 
clouds. 

We focus our study initially on the p Oph region. Our 
goal is to construct a radiative transfer model to estimate the 
temperature of the cores in this region taking into account 
the 3D geometry of the region and the role of HD 147879 in 
externally heating the clumps. We will discuss the applica- 
bility of this study to other regions of star formation where 
low-mass stars form, and to more energetic regions of star 
formation, such as the Orion nebula, by discussing the role 
of a stronger external radiation field. 

The structure of the paper is as follows. In Section 3 
we discuss the radiative transfer method we use and the 
constituents of the model (i.e. radiation sources, model ge- 
ometry, dust properties, clump density profile). In Section 4 
we describe the results of the radiative transfer modelling 
with respect to the temperature profile and the SEDs of the 
clumps in p Oph, and in Section 5 we discuss the implication 
of the estimated dust temperature for the core mass function 
of p Oph. In Section 6, we discuss the dust temperatures of 
cores in other star forming regions, and finally in Section 7 
we summarize our results. 



3 THE MODEL 

The big clumps of p Oph are represented by Plummer-like 
density profiles, i.e. we invoke a spherical geometry, with the 
density being approximately flat in the centres of clumps 
and dropping as in the envelopes. We exclude the Oph- 
Bl clump, which is quite flattened, and the Oph-E clump, 
which appears to be part of Oph-C. We shall assume that the 
cores in the clumps are superimposed on this density profile. 
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n=n„/[H-(r/rof] 

Dust: MRN standard mixture 
+grain coagulation 
+thin ice mantles 



Ambient cloud 



HD 147889 




p= constant 

Dust: MRN standard mixture 



Observer (6-120") 



Figure 1. Schematic representation of the p Oph A model. The p 
Oph A clump is surrounded by a virtual ambient cloud that mod- 
ulates the incident radiation that heats the clump. The external 
heating is provided both by the interstellar radiation field and by 
HD147889, a nearby B2V luminous star. The general geometry 
of the models for the other p Oph clumps are similar to this one. 
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Figure 2. Dust opacity for the clumps (Ossenkopf & Henning 
1994) and the virtual ambient cloud (Draine 2003). The solid lines 
correspond to the opacity due to absorption, and the dashed lines 
to the opacity due to scattering. 



without actually taking them individually into account in 
this model. However, in Section 14.71 we shall estimate their 
effect on the dust temperatures. We also assume that these 
clumps are heated (i) by the interstellar radiation field and 
(ii) by radiation from the HD 14789 B2V star (see FiglU. We 
will adopt the 3D model constructed by Liseau et al. (1999). 
Hence, in our model these clumps define the optical depth 
that the external radiation has to penetrate to reach the 
individual cores embedded inside the clumps. In the next 
subsections we describe in detail the model. 

3.1 Density profiles 

Each of the p Oph clumps is represented by a density profile, 
'^('^) = no , (2) 

where no is the density at the centre of the clump, and ro 
is the extent of the region in which the density is approxi- 
mately uniform. This profile is approximately fiat in the cen- 
tral region and falls of as in the envelope, similarly to the 
more commonly used Bonnor-Ebert sphere density profile 
(Ebert 1955; Bonnor 1956). Additionally, despite being ad 
hoc, it fits well the observed profiles of prestellar cores and 
it predicts hfetimes, accretion rates, collapse velocity fields, 
SEDs and isophotal maps which agree well with observation, 
using a minimum number of free parameters (Whitworth & 
Ward- Thompson 2001). The adopted values for no, ro and 
the extent 7?ciump of each clump are listed in Table [T] These 
values are consistent with observations of these clumps (e.g. 
Motte et al. 1998; Andre et al. 1993). We assume that the 
cores in the clumps are superimposed on this density profile 
(we re-iterate that at this stage we ignore the sub-structure 
of the clump). 

Each clump is surrounded by a virtual ambient cloud, 
which has a uniform density. The role of these ambient 



clouds is to modify the ambient radiation field that heats 
each clump externally. The optical extinction through the 
ambient cloud is chosen so that the peak of computed SED of 
each clump corresponds to a dust temperature that matches 
the dust temperature of each clump assumed by previous au- 
thors. The optical extinctions of the assumed ambient clouds 
around each one of the p Oph clumps are also listed in Ta- 
ble [1] 



3.2 Dust properties 

We use two kinds of dust; one for the clumps, and one for the 
ambient cloud (Fig. [2| . The dust grains in dense clumps are 
expected to coagulate and accrete ice mantles, so we use the 
Ossenkopf & Henning (1994) opacities for MRN dust that 
has coagulated and accreted thin ice mantles for a period 
of 10''' years at a density lO'^ cm~^. The ambient cloud has 
relatively low density and it can be considered as a part 
of the interstellar medium. Hence, for this layer a standard 
MRN dust opacity is used (Draine 2003). 



3.3 Radiation sources & geometry of the model 

Each clump is heated (i) by the interstellar radiation 
field, attenuated through the ambient cloud, and (ii) by 
HD147889, a B2V nearby star which is located at distance 
~ 0.5 — 1 pc from the p Oph region (e.g. Liseau et al. 1999). 
In this model we will ignore radiation from deeply embedded 
protostars in these clumps (e.g. VLA 1623, IRS43). Accord- 
ing to the Stamatellos et al. (2005) model the presence of an 
embedded protostar heats only the region around it, hence 
the results presented here are not greatly affected by the 
presence of young protostars in or around the clumps. In 
Section 14.61 we shall discuss the effects of such embedded 
protostars. 
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Table 1. (I) Assumed parameters for each clump (based on observations). (II) Results of the radiative transfer modelling. 



Parameter p Oph A p Oph B-2 p Oph C p Oph D p Oph F 





-Rclump ^ (AU) 


2 X 10* 


2.9 X 10* 


1.7 X 10* 


1.2 X 10* 


1.2 X 10* 




no ^ (cm~'^) 


16 X lO** 


1.2 X 10"^ 


1.6 X lO*' 


5.2 X 10^ 


3.3 X 10^ 




r-o J(AU) 


10=* 


5.4 X 10^ 


3.4 X 10» 


4 X 103 


7 X 10^ 


(I) 


4 cloud d 

— 


0.19 


0.35 


2.8 


0.45 


0.17 




d*l(pc) 
M /. (Mq) 


1.1 


1.4 


1.3 


1.8 


1.6 




24 


61 


18 


47 


14 




^clump g 


1090 


400 


560 


120 


130 




rcentroil(K) 


6 


7 


7 


7 


8 


(11) 


rprnin i (j^) 
edge ™ 


14 


13 


12 


14 


14 




edge " ■' 


27 


22 


14 


20 


24 




Apeak 


125 


200 


270 


200 


160 



Clump extent. 
^ Clump central number density. 

Flattening radius of each clump. 

Visual extinction through the ambient cloud. 
^ Distance of each clump from HD 147889 

f Clump mass (calculated from the density profile assumed and the extent of the clump). 

9 Visual extinction to the centre of each clump (calculated from the density profile assumed and the extent of the clump). 
^ Temperature at the centre of the clump. 

* Temperature at the edge of the clump facing away from the H147989 direction. 
J Temperature at the edge of the clump facing towards H147989. 
^ Wavelength of the peak of the computed SED 



For the stellar radiation we adopt the parameters com- 
puted by Liseau et al. (1999) {Tf = 22, 000 K, ii* = 5 Rq, 
Li, — 5, 300 Lq). We shall assume that this star emits as a 
blackbody having temperature 22,000 K. Due to the pres- 
ence of this star in the vicinity of p Oph, the external heating 
of the clumps is highly anisotropic. 

For the external interstellar radiation field we adopt 
a revised version of the Black (1994) interstellar radiation 
field (BISRF). The BISRF consists of radiation from giant 
stars and dwarfs, thermal emission from dust grains, cosmic 
background radiation, and mid-infrared emission from tran- 
siently heated small PAH grains (Andre et al. 2002). This 
radiation field is modulated by the ambient cloud around 
each clump, hence the incident radiation field on each of 
the clumps is enhanced at FIR and longer wavelengths, and 
attenuated a shorter wavelengths. 

The geometry of the p Oph region is taken from the 3D 
model constructed by Liseau et al. (1999) based on far in- 
frared spectrophotometric observations with the ISO-LWS. 
Using this model and the 1.3mm mosaic image of Motte et 
al. (1998), we estimate the distance of HD 147889 from each 
of the p Oph clumps (see Table[T|| . These are rough estimates 
but they do not greatly affect the results of the model. 

3.4 Monte Carlo radiative transfer 

The radiative transfer calculations are performed using 
PHAETHON, a 3D Monte Carlo radiative transfer code 
developed by Stamatellos & Whitworth (2003). The code 
uses a large number of monochromatic luminosity packets 



to represent the radiation sources in the system. The lumi- 
nosity packets are injected into the cloud and interact (are 
absorbed, reemitted, scattered) stochastically with it. If an 
L-packet is absorbed its energy is added to the local re- 
gion and raises the local temperature. To ensure radiative 
equilibrium the L-packet is re-emitted immediately with a 
new frequency chosen from the difference between the local 
cell emissivity before and after the absorption of the packet 
(Bjorkman & Wood 2001; Baes et al. 2005). 

The model for each p Oph clump is essentially 2- 
dimensional, hence the code used here is adapted and op- 
timised for the study of systems with azimuthal symmetry. 
Each clump is divided into a number of cells by spherical and 
conical surfaces. The spherical surfaces are evenly spaced in 
radius, and there are typically 100 of them. The conical sur- 
faces are evenly spaced in polar angle, and there are typically 
40 of them. Hence, each clump is divided into ~ 4000 cells. 
The number of cells used is chosen so that the density and 
temperature differences between adjacent cells are small. 

The L-packets representing the ambient radiation field 
(typically a few 10^" packets) are injected from the outside 
of the clump with injection points and injection directions 
chosen to mimic an isotropic radiation field incident on the 
ambient cloud around each clump (Stamatellos et al. 2004). 

The L-packets representing the stellar radiation from 
HD147889 (typically a few 10^'' packets are used) are emit- 
ted from the star with random direction and only a fraction 
of them (~ 1%) heats the clumps. However, due to the large 
luminosity of HD147889, this small percentage of stellar pho- 
tons dominates over the background radiation field. 
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4 DUST TEMPERATURES AND SEDS OF 

THE CLUMPS IN THE p OPH MAIN CLOUD 

The optical depth of the virtual ambient cloud around each 
clump is varied so that to produce an SED for each clump 
that peaks at a wavelength corresponding to a dust temper- 
ature that matches the dust temperature assumed by pre- 
vious authors. In the next subsections we present the dust 
temperature profiles and the SEDs of each clump as calcu- 
lated by this model (see Figs. I3I7|| . 

4.1 p Oph A 

Oph A is, according to the Liseau et al. (1999) 3D model, 
the closest clump to HD147889. The radiation from the star 
dominates over the interstellar radiation field and heats the 
clump to temperatures that range from 6 — 7 K in the cen- 
tre of the clump, to 13 — 27 K at the clump edge (Fig. |3]). 
The clump hemisphere closer to the star has higher tem- 
peratures (up to 27 K) than the other hemisphere, which 
has temperature not higher than 15 K. In total, 80% of the 
clump volume is colder than 15 K. The ambient cloud has 
dust temperatures 30-40 K on the side closer to the star, 
which are consistent with previous temperature estimates 
based on the PDR model of Liseau et al. (1999). 

The SED of the model of p Oph A (Fig. |31d) peaks 
at around 125 ^.m, which is suggestive of a temperatur^fl 
around ~ 20 K. This value is consistent with what other 
studies have indicated (e.g. Andre et al. 1993) about the 
temperature derived from the SED. However, the peak of 
the SED peaks at this wavelength due to the contribution 
from the outer hotter parts of the clump, whilst most of the 
clump is colder. The region of the clump where prestellar 
cores are observed, i.e. at visual extinctions Ay > 7 mag 
(according to Johnstone et al. 2004), is colder than 11 K. 

4.2 p Oph B-2 

The temperature profile of the p Oph B-2 clump (Fig. |4]) is 
similar to that of Oph A. This clump is colder than Oph 
A, and its temperature drops from 13 — 23 K at the edge 
of the cloud to ~ 7 K at its centre. The SED distribution 
of the clump peaks at around 195 \im indicating a tempera- 
ture of 13 K, which is consistent with previous observations. 
As in Oph A, the SED peak is characteristic of the outer 
warmer layers, whereas most of the clump is colder than 
13 K. Again assuming that cores are observed at visual ex- 
tinctions Av > 7 mag, we find that their temperatures 
should be below 10 K. 

4.3 p Oph C 

The estimated temperature of this clump from previous ob- 
servations is around 10 K (e.g. Motte et al. 1998; Andre et 
al. 1993). This clump is colder than the previous two, sug- 
gesting that it is shielded from HD147879, possibly by parts 
of the other p Oph clumps. The temperature of this clump 

We note that the peak of SED corresponds to the peak of 
A kxBx, which for kx oc i.e. the Ossenkopf & Henning 

(1994) opacity, translates to App^j^T 2, 500 K [am. 
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Figure 3. (a) Dust temperature profile of the Oph-A clump ver- 
sus the distance from its centre, and versus visual extinction from 
the surface of the clump. We plot the radial temperature profile 
from 9 = 0° to 180° to every 9° , where 9 = 0° corresponds to the 
direction towards the B2V star. The red thick line corresponds to 
the direction perpendicular to the clump-star direction (red thick 
strip in Fig. [T}. (b) Simulated SED of the clump. 

drops from 11-13 K at its edge of to 6-7 K at its centre 
(Fig. [5^). The SED of the system peaks at around 270 yun 
(Fig. [S]d), which corresponds to a temperature of 9 K. As- 
suming that cores exist at visual extinctions Ay > 7 mag 
we find that their temperatures should be below 10 K. 

4.4 p Oph D 

The temperature profile and the SED of p Oph D are pre- 
sented in Fig. |6l The SED peaks at around 205 \im, corre- 
sponding to a temperature of 12 K, as assumed by Motte et 
al. (1998). However, most of the clump is colder than 12 K. 
Assuming that cores exist at visual extinctionSj4v > 7 mag 
we find that their temperature should be below 10 K. 

4.5 p Oph F 

The temperature profile and the SED of p Oph F are pre- 
sented in Fig. [T] The SED peaks at around 160^m, corre- 
sponding to a temperature of 16 K. However, most of the 
clump is colder than 16 K, with the region at visual extinc- 
tions Ay > 7 mag having temperature below 10 K. 

4.6 Dust heating by young protostars 

We have so far ignored heating from nearby young proto- 
stars. Motte et al. (1998) identify a number of Class 0, Class 
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Figure 4. Same as in Fig. |3] but for the Oph-B2 clump. 
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Figure 5. Same as in Fig.[3]but for the Oph-C clump 
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Figure 6. Same as in Fig. |3]but for the Oph-D clump. 
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Figure 7. Same as in Fig. |3] but for the Oph-F clump. 



The dust temperatures of prestellar cores 7 



10-' e 




Figure 8. SED of a young Class object as calculated by Sta- 
matellos et al. (2005), assuming a distance of 160 pc. The bolo- 
metric luminosity of this source is 5 Lq. 
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Figure 9. (a) Dust temperature profile of the clump, in different 
directions, due to heating from a Class object at the edge of the 
clump (b) Dust temperature profile, in different directions, due to 
heating from a Class object and the interstellar radiation field. 
The thick (red) line corresponds to the direction perpendicular 
to the clump-protostar direction. The dashed (blue) line corre- 
sponds to the temperature when only heating from the ISRF is 
considered. The role ISRF is dominant. 



I and Class II objects either near the clumps or embedded 
in the clumps. The presence of these luminosity sources is 
expected to increase the clump temperature. In order to 
quantify this increase we use the p Oph A clump as a study 
case, and consider heating (i) from a Class object located 
at the edge of the clump, (ii) from a Class I object also lo- 
cated at the edge of the clump, and (iii) from a Class I object 
embedded in the clump, 0.05 pc from its centre. The SEDs 
of the Class and the Class I objects are taken from the 
simulations of Stamatellos et al. (2005) (Figs. |S1 [Tg. The 
bolometric luminosity (i.e. the integrated luminosity over 
the entire spectrum) of the Class source is 5 Lq, and that 
of the Class I source is 10 Lq . These are typical luminosities 
of such objects in the region (e.g. Wilking et al. 1989). 

The results of the radiative transfer simulations are pre- 
sented in Figs.l^ im and ll2l On the top graph of each figure 
we present the dust temperature taking into account only 
the heating from the nearby or the embedded source. The 
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Figure 10. SED of a young Class I object as calculated by Sta- 
matellos et al. (2005), assuming a distance of 160 pc. The bolo- 
metric luminosity of this source is 9.3 Lq. 
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Figure 11. Same as in Fig. [9] but for heating from a Class I 
object at the edge of the clump. The ISRF is again dominant. 



clump is hotter closer to the assumed source, as expected, 
with the temperature ranging from 4 to 15 K in different 
directions inside the clump. On the bottom graphs of each 
figure we present the dust temperature also taking into ac- 
count the heating of the ISRF, which has a bolometric lumi- 
nosity of 5 Lq, i.e. comparable with the luminosities of the 
nearby or embedded sources. As seen in the figures the role 
of the ISRF is dominant; the presence of the nearby Class 
object does not affect the dust temperature at all, and 
the presence of a nearby Class I object just heats the outer 
layers of the clump that are closer to it by only ~ 2 K. The 
presence of an embedded Class I object is more noticeable. 
It increases the dust temperature in the clump by ~ 1 — 2 K, 
and very close to it by ~ 5 K. However, we note that (i) the 
temperature of most of the clump is still below 10 K, and 
(ii) the increase of temperature is not greater than the in- 
crease of temperature when the radiation from HD 147889 is 
considered (cL Fig.|3]D). This is because (i) the luminosity of 
the protostar is comparable to the luminosity of the ISRF, 
hence smaller than the luminosity of HD 147889 reaching the 
clump, and (ii) most of the protostar's radiation is emitted 
at long wavelengths, where the clump is optically thin. 

We conclude that the effect of heating due to the pres- 
ence of nearby or embedded young protostars is not greater 
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Figure 12. Same as in Fig. |9] but for heating from a Class I 
object located inside the clump (0.05 pc away from the centre of 
the clump). The heating provided by the ISRF is still dominant 
apart from the region close to the embedded protostar. 



than the effect of the ISRF, and hence secondary to the effect 
of HD147889. Thus, nearby or embedded young protostars 
do not greatly affect the dust temperatures in the regions 
deep inside the clumps where prestellar cores are observed, 
unless there is a large number of young protostars within 
each clump, which is not the case for the p Oph main cloud. 

4.7 The effect of a 3D clumpy structure 

In the radiative transfer simulations presented so far we use 
spherically symmetric models to describe the big clumps in 
the p Oph main cloud. In reality these clumps are more 
structured, as they contain regions where the density is ei- 
ther higher or lower than the assumed spherically symmet- 
ric model. Indebetouw et al. (2006) examined the effect the 
cloud dumpiness in the case of clouds that contain high- 
mass stars and they found that the dust temperature is af- 
fected significantly due to the dumpiness of the medium. 

To examine the effect of the 3D clumpy structure of 
the regions of the p Oph main cloud, we generate a 3D p 
Oph-A-like dump as follows. We start off with a spherically 
symmetric clump with the parameters adopted for p Oph-A 
(see Table [l]). We then impose a turbulent velocity field (e.g. 
Goodwin et al. 2004), and use DRAGON, an SPH code, to 
follow the evolution of the cloud. Due to the effect of tur- 
bulence the cloud acquires a clumpy structure. The cloud is 
evolved until cores are formed (Figs. [131 114p . We then per- 
form a radiative transfer simulation on the clumpy structure 
using the method of Stamatellos & Whitworth (2005a). As 
in the previous sections we consider heating both by the 
ISRF and HD147889. The calculated temperature profile is 
presented in Fig. 1141 The dust temperature is similar to the 
temperature calculated in the spherically symmetric case, 
despite the fact that the density in some regions of the cloud 
is up to an order of magnitude different. This is because the 
heating of the dense parts of the cloud is mainly due to long 
wavelength radiation that propagates into the cloud without 
"seeing" the clumpy structure; short wavelength radiation is 




I 



Figure 13. Logarithmic density (in g cm~^) on the x — y plane 
of a simulated Oph A-like clump. The density ranges from 10~^* 
to lO^^'' g cm^''. X and y are given in cm. 
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Figure 14. Density and temperature profiles (dots) of the sim- 
ulated clump of Fig. 1131 Each dot corresponds to the tempera- 
ture of each radiative transfer cell constructed from the output of 
the SPH simulation. The density and temperature profiles of the 
spherically symmetric model are also plotted (lines). The upper 
green line corresponds to the direction towards HD 147889 and the 
lower green line to the opposite direction. The dust temperature 
in the case of a 3D clumpy structure is similar to the temperature 
calculated using the spherical symmetric model. 



absorbed in the outer layers of the cloud and re-emitted at 
longer wavelengths. 

We conclude that the dumpiness of the cloud does not 
significantly affect the dust temperature at the inner, dense 
regions of the cloud. 
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4.8 The dust temperature of the cores in the p 
Oph main cloud 

According to the model presented here, the dust tempera- 
tures of the cores in p Oph main cloud are lower than previ- 
ously thought. Since cores are observed at visual extinctions 
> 7 mag, the temperatures of the cores in Oph-A are most 
probably below 11 K, and in Oph-B2, Oph-C, Oph-D and 
Oph-F below 10 K. These temperatures correspond to the 
temperature of the clumps where the cores are embedded, 
at visual extinctions Av = 7 mag. 

Considering the fact that most of the cores are observed 
at visual extinctions 12 — 20 mag (Johnstone et al. 2004), 
their actual temperatures may be even lower by ~ 2 K. 
Hence, the presence of the luminous HD 147889 in the vicin- 
ity of the clumps, affects only the outer regions of the 
clumps, where prestellar cores have not been observed. 
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Figure 15. The core mass function of the p Oph main cloud 
based on the core masses calculated by Motte et al. (1998). The 
error bars correspond to \/N counting statistics 



5 THE CORE MASS FUNCTION OF p OPH 

Motte et al. (1998) calculated the core mass function of the 
p Oph main cloud by assuming a temperature of 20 K for 
Oph-A, 12 K for Oph-B, Oph-C, and Oph-D, and 15 K for 
Oph-E and Oph-F. Johnstone et al. (2000) have assumed 
similar temperatures to Motte et al. (1993) or even higher 
temperatures. These temperatures are higher than the ones 
predicted by our model, by 2 — 9 K. Hence, considering the 
fact that core masses are calculated using Eq. [T] we suggest 
that both these authors have underestimated the masses of 
the prestellar cores in this region. At these temperatures the 
Rayleigh- Jeans relation is not a good approximation and the 
Planck function must be used; hence even overestimating 
temperatures by a few K leads to underestimating masses 
by a factor of 2 — 3. 

We revise the core mass function of Motte et al. (1998) 
assuming core temperatures derived from our model (Sec- 
tion |4]8]). We shall assume the same dust opacity at 1.3 mm, 
i.e. K = 0.005 cm^ g~^ (Preibisch et al. 1993; Andre et al. 
1996). We then calculate the core the core masses from the 
observed fluxes using Eq [T] We note again that tempera- 
tures we assume correspond to cores at visual extinctions of 
~ 7 mag. Johnstone et al. (2004) suggest that most of the 
cores in this region are embedded at 12 — 21 mag or at even 
higher extinctions. Hence the temperatures of many of the 
cores are expected to be even lower (by ~ 2 K). This means 
that the temperatures used are upper limits to the actual 
temperatures, and consequently that the computed masses 
are lower limits to the actual core masses. 

In Fig.[TS]we present the Motte et al. (1993) CMF, and 
in Fig.[T6]the revised CMF derived using core temperatures 
from the model presented in this paper. The CMF has moved 
to higher masses but its overall shape has not changed signif- 
icantly. The mass at which the core mass spectrum steepens 
from a slope a ~ 1.5 to a slope a ~ 2.5 is less clear in 
this case but it seems that it has moved from ~ 0.5 Mq to 
~ 1 Mq. Contrary to the CMF in other star forming re- 
gions (e.g. in Orion; Nutter et al. 2007) the CMF does not 
show a turnover down to the completeness limit (which is 
^ 0.2 Mq, using the new temperatures). However, the CMF 
may flatten at around ~ 0.4 Mq. 
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Figure 16. The CMF of the p Oph main cloud based on the 
core masses calculated using temperatures from our model. The 
CMF is constructed so that there is approximately equal number 
of cores per bin. The CMF has moved to larger masses but its 
overall shape has not changed significantly. 



6 DUST TEMPERATURES OF CORES IN 
OTHER STAR FORMING REGIONS 

The dust temperatures of prestellar cores are determined 
by the ambient radiation field heating the core (since by 
definition there are no radiation sources inside the core). 
This radiation field is determined by environmental factors. 
It depends (i) on how deeply the core is embedded in its 
ambient cloud, and (ii) on the possible presence of nearby 
radiation sources, and their relative position/distance with 
respect to the core. Hence, the external radiation field that 
heats each core is different. 

Previous studies of prestellar cores and young proto- 
stars (e.g. Evans et al. 2001; Young et al. 2003; Jorgensen 
et al. 2006) have acknowledged this fact and have used a 
scaled version of the standard ISRF that is either enhanced 
at all wavelengths or selectively at UV and FIR. This sim- 
ple approach has a free parameter, the ISRF scaling factor, 
which is varied arbitrarily to fit the observations but it is not 
connected directly to the molecular cloud in which the core 
is embedded or the transport of radiation inside the cloud. 
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It also does not account for the fact that the radiation field 
incident on an embedded core is not isotropic. 

Here, we follow the approach of Stamatellos et al. (2004) 
where the radiation field incident on the core is a direct 
result of the presence (i) of the ambient cloud that sur- 
rounds the core, (ii) of nearby embedded young protostars, 
(iii) of nearby luminous stars. The ambient cloud of each 
core/clump attenuates the ISRF, as it acts like a shield to 
UV, visual and NIR interstellar radiation, absorbing and 
re-emitting it in the FIR. The nearby young protostars 
and stars also enhance the radiation field incident on core, 
mainly at long wavelengths (> 50 \im) due to the thermal 
emission from the outer regions of the heated ambient cloud 
(Mathis et al. 1983), as short wavelength radiation cannot 
penetrate deep into the cloud. 



6.1 Dust temperatures of the prestellar cores in 
the p Oph molecular cloud 

In the previous sections we discussed in detail the region of 
the p Oph main cloud, where the external heating is dom- 
inated by a nearby B2V star. Due to the presence of this 
star and embedded protostars in the p Oph main cloud the 
external radiation field is stronger than the standard ISRF 
by a factor of ~ 10 (the bolometric luminosity of the ISRF 
heating the core is 5 Lq whereas the bolometric luminos- 
ity of the star's radiation reaching and heating each clump is 
~ 50 Lq). This is consistent with the observations of Liseau 
et al. (1999) which suggest that the external radiation field 
incident on p Oph ~ 10 — 100 times stronger than the stan- 
dard ISRF. 

Using a detailed model presented in the previous sec- 
tions, we estimate that the temperatures of the cores in the 
p Oph main cloud are probably below 10-11 K, i.e. lower 
than previous estimates. The dust temperatures of the cores 
in L1689, another region of the Ophiuchus molecular cloud 
which is totally starless (e.g. Nutter et al. 2006), are also 
probably below 10-11 K. 



6.2 Dust temperatures of the prestellar cores in 
the Taurus molecular cloud 

In Taurus there are no nearby luminous stars, hence the ra- 
diation field heating the cores in this region is enhanced only 
due to the presence of young protostars and stars in its vicin- 
ity. This radiation field is not expected to be stronger than 
that in p Oph. Thus, the dust temperatures estimated by 
our model for the p Oph main cloud are upper limits to the 
temperatures of the cores in Taurus. Indeed, observations of 
TMC-IC, a prestellar core in in Taurus (Schnee & Good- 
man 2005; Schnee et al. 2006), suggest dust temperatures 
that drop from 11 K at the edge of the core {Ay « 10) to 
5 K at its centre {Ay « 80), which are consistent with heat- 
ing from an ISRF that is weaker than the standard ISRF. 

In Fig. [17] we present a model for TMC-IC. Based 
on the estimates from Schnee et al. (2006) we assume 
the density profile defined in Eq. [S] with r,, = 0.02 pc, 
no = 5.5 X 10^ cm~^ and Ticoro = 15000 AU. We further 
assume a virtual ambient cloud around the core of visual 
extinction Ay = 3.5 . This virtual cloud modifies the radia- 
tion field incident on the core. The calculated temperature 
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Figure 17. Assumed density profile and calculated temperature 
profile of TMC-IC (solid lines) in Taurus. The squares corre- 
spond to estimates based on observations at 450 yim, 850 yim, 
and 1200 lam (Schnee et al. 2007). 



(Fig. 1171 bottom) is very similar to the temperature esti- 
mated by Schnee et al. (2006). The calculated temperature 
at the centre of the core is higher by ~ 1 K than the Schnee 
et al. (2006) estimate, indicating that the radiation field 
heating the core is weaker than the standard ISRF at long 
wavelengths (or alternatively that the core centre is denser) . 

We conclude that the dust temperatures of the prestel- 
lar cores in Taurus are probably below 10 K. 



6.3 Dust temperatures of the prestellar cores in 
the Orion molecular cloud 

In Orion there are many radiation sources that contribute 
to the heating of the prestellar cores in region. Jorgensen 
et al. (2006) estimates that the external radiation field in 
the Orion region is up to ^ 10^ stronger than the standard 
interstellar radiation field. 

In order to investigate the effect of such an enhanced ra- 
diation field on the dust temperatures in prestellar cores in 
Orion, we use p Oph A as a representative clump, and con- 
sider heating from a B2V star (similar to HD 147889) with 
different temperatures and at different distances from the 
centre of the clump. Hence, we examine the effect of heating 
from an incident radiation field having bolometric luminos- 
ity Lh 10, 20, 50, 120, 500 and 10^ times the bolometric 
luminosity of the standard interstellar radiation field (Black 
1994; Andre et al. 2002). 

In Fig. [18] we present the dust temperature versus vi- 
sual extinction calculated using incident radiation fields of 
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Figure 18. Dust temperature of a p Oph A-like clump against the visual extinction from the edge of the clump. The clump heating is 
due to the ISRF and a B2V star having different distances from the centre for the clump and different temperatures, so as to change the 
bolometric luminosity incident on the clump (as marked on each graph). The temperature of most of the clump is below 10 K unless the 
heating from the star is very strong (iJ 120 X Lisrf). 



different strengths. For radiation fields with bolometric lu- 
minosities up 50 X LisRF the dust temperature profile is 
similar to the ones calculated for p Oph; most of the dust 
at optical visual extinctions > 7 mag is colder than 10 K. 
For a radiation field 120 x Lisrf, half of the clump (i.e. the 
hemisphere away from the luminosity source) is colder than 
12 K. For a higher external radiation, (500 — 10^) x Lisrf, 
the dust temperature in the clump varies considerably with 
the position in the clump; from 14 to 45 K in the case of 
500 X Lisrf, and from 15 to 53 K for the case of 10^ x Lisrf. 
Hence, the dust temperature of a core depends on the core 
position in the clump, i.e. its relative position with respect to 
the external radiation source and on how deeply embedded 
the source is in its parent cloud. 



We conclude that assuming that (i) the incident radi- 
ation field on the clumps of Orion is (500 — 10"^) x Lisrf, 
and (ii) that cores exist at visual extinctions > 7 mag, then 
typical dust temperatures of the prestellar cores in Orion 
are from 20 to 30 K, which is consistent with previous as- 
sumptions/estimates (Launhardt et al. 1996; Mitchell et al. 
2001; Motte et al. 2001; Johnstone et al. 2006; Nutter et al. 
2007). 



7 CONCLUSIONS 

We have used a 3D geometry for the p Oph main cloud 
region to construct a radiative transfer model for this re- 
gion, taking into account external heating by (i) the inter- 
stellar radiation field, and (ii) HD147889, a nearby B2V star. 
HD147889 dominates the heating of the p Oph main cloud 
clumps. We estimate that the dust temperatures at visual 
extinctions > 7 mag are below ~10-11 K. These are smaller 
than was previously assumed. As a result we find that the 
core masses calculated from mm observations are underesti- 
mated by a factor of 2-3. This affects the core mass function 
of the p Oph main cloud. We computed a revised CMF for 
this region using the dust temperatures calculated in this 
paper. The CMF has moved to higher masses but its shape 
has not changed significantly. The mass at which the core 
mass spectrum steepens from a slope a ~ 1.5 to a slope 
a ~ 2.5 is less clear but it appears that it has moved from 
~ 0.5 Mq to ~ 1 Mq. This is still below the mass where 
this steepening occurs in Orion. Contrary to the CMF in 
Orion the CMF of the p Oph main cloud does not show a 
turnover at low masses. However, it may flatten at around 
~ 0.4 Mq. 

We have generalized our study to estimate the dust tem- 
peratures in prestellar cores in other star forming regions. 
In the Taurus molecular cloud the ambient radiation field is 
weaker than that in Ophiuchus, hence the dust temperatures 
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of the cores in this region are similar or smaller than the ones 
calculated for p Oph. We estimate that the dust tempera- 
tures at visual extinctions > 7 mag are below ~10 K. In 
Orion the ambient radiation field is estimated to be up to 
10'^ times stronger than the standard interstellar radiation 
field. Using a simple model to account for this enhanced ra- 
diation field, we find that the typical dust temperatures of 
the prestellar cores in this region are from 20 to 30 K. 
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